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Abstract 


tow^S»,kTw yd , r0 f , ™ , " iC (MHD) ™ 8,ne,,led pro8ul5io " reveak • operating regime a. 

, . 1 s - Substanti al mass savings are realized over conventional chemical, solar, and electrical 

p opulsion concepts when MHD augmentation is used to obtain optimal I sp . However trip times for the 

Zl-t° nS , e T Ve o tCS ° f P ° Wer Plant specific N»he and accelerator efficL™ may be 
Zm b V a k x'lfm Quasi-one-dimensional calculations show that a solar or nuclear thernal system 
gmented by MHD can provide competitive performance while utilizing a diverse range of propellants 

lt.K ng r ’ Wh ‘, Ch ,S aVaUable fr ° m the SpaCe Shuttle ’ the droids, and vafious moons Td 

as enahl h hm ^ ^ ^ USe of in ' sltu P ro P e,Iants wiI1 reduce costs of space operations as well 

as enable human exploration of our Solar System. 

The following conclusions can be drawn from the results of the mission trade study: 

' raaSta ° r maSS * , °' V * ab ° Ve MHD ““ 

Mass saving of over 50% can be realized for unique combination of soIar/MHD systems 

Trip times for systems utilizing cuirent power supply technology may be prohibitively long. 

p ° ret,cal P red 'ctions of MHD performance for in space propulsion systems show that improved 
efficiencies can reduce trip times to acceptable levels. 

Long trip times indicative of low thrust systems can be shortened by an increase in the MHD 
accelerator efficiency or a decrease in the specific mass of the power supply and power processing 

- As for all propulsion concepts, missions with larger Av’s benefit more from the increased specific 
impulse resulting from MHD augmentation. p c 


a quas ‘? nedimensional analysis, the required operating conditions for a MHD accelerator to reach 
acceptable efficiencies are outlined. This analysis shows that substantial 
desirable. 


non-equilibrium ionization is 


Introduction 


MHD augmentation has been discussed, analyzed, and experimentally evaluated for over 50 years. In fact 
artmann carried out the first investigation of MHD channel flow in the 1930s, examining MHD effects in 
a steady incompressible conducting fluid. Beginning in the late 1950s, MHD was looked at to provide 
electrical power, hypersonic testing facilities, drag reduction, air flow manipulation, and augmentation to 
propu sion systems. It is estimated that over 1000 man-years have been spent analyzing MHD effects since 


In recent years there has been a renewed interest in MHD for use in hypersonic testing facilities. Igor V. 
Adamovich and J. William Rich of Ohio State University working with Gordon L Nelson of MSE Inc. 
completed a feasibility study on MHD acceleration in seeded and unseeded airflow. In the unseeded 

fnlff W t’ an ! eCtr °? b r m 1 s I ustamed lonization - At hi gb pressures electron beam ionization was shown to be 
ineffective due to high collision frequencies. However, at lower pressures non-equilibrium ionization was 
effective although test section pressures would be too low to allow its use in hypersonic wind tunnels. 

wm extemal| y sustained ionization will not work for wind tunnel applications it shows promise for use in 
MHD augmented upperstage propulsion systems. H 

Several recent reports and papers have been published on the feasibility of using MHD augmentation for 
upperstage propulsion systems. Schulz and Chapman have shown that a MHD augmented chemical 
dmister can achieve specific impulses up to 4000 m/sec. Bagher M. Tabibi et. al. have explored the use of 
MPD augmentation of solar propulsion systems. To date, no system level study has been completed to 
determine the optimal operating regime for MHD augmented systems. 
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nntiliT eI r CtnC Pr T' Si0n SySt6m ’ MHD au 8 me nted systems have an optimal specific impulse Tire 

nro^Zt c CC1 hVT 15 3 re f Uk ° f the maSS 0f the P° wer su PP'y increasing faster than the mass of file 
p p ant sa ved by the increased specific impulse. Figure 1 shows the relationship between the initial mass 

of a spacecraft, the power supply mass, and the mass of the propellant. 




Specific ImpuEie (see) 


In itial M ass (M T) 

Propellant Mass (MT) 
Power Plant M ass (M T) 


Figure 1 : The relationship between the initial mass, power plant mass, and propellant mass. 

The kinetic power of the exhaust jet is given by 

P jet(T.I S p) ^ gT I sp (1) 

ZX T n fVh h . e ^! ' ' 0Ck “' g ; he »f SnmV at the eanh's surface, and I„ is the specific 

impulse of the rocket. The corresponding mass flow rate is then given by P P 

m 'XT'! 2 Pjet( T ’ Is p) 

m dot( T » I sp) - — — 1 ,3 (2) 

( 8 sp) 

Keeping the mass flow rate constant, the jet power of the MHD system for a given bus power and 
efficiency of conversion of bus power to kinetic energy in the exhaust (p) is determined by 

P jet_MHD( T >Isp.Pbus>n) :=p jet(T,I sp ) + P b us1 (3) 

Then the specific impulse of the MHD augmented system is simply 


r t t n \ . 1 2 P jet MHDf r > r sp> p bus> T l') 

>sp_ M m>T,i sp .P bus , n ) ~.J _J 2 


dot > sp) 
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Likewise the thrust of the MHD augmented system is obtained by 

T MHD T^sp-Pbus’lj ■ -m dot ; i T,I sp) g' I sp_MHD l T’ I sp’ P bus> r l | (5) 


m 


elect( a > p bus) :_aP bus 


( 6 ) 


Where a is the specific mass of the power system. In this study the mass of the power-processing unit has 

“ SP ' I0 ° nUClM, P”™ sy *' m “ available to deliver the 
tl! P m3SS ° f thlS SyStem ,s 12 k § /kW - ^ s P ecific mass of the power processing 

d conditioning system was set at 7 kg/kW giving a = 19 kg/kW. The affect of decreasing the specific 
mass of the power supply (a) will also be looked at. 8 P 

The mass of the solar thermal system can be broken up into two parts the mass of the mirrnrc 

rocket engine - For a ,ightweighi continuous thrust system the 


m STR^) := 0-!25T 


1. 15 


( 7 ) 


tl ) rust ; T ’ is m Newtons. Marshall Space Flight Center has also used the following scaling for 
solar thermal systems, which utilize a store and bum technique. 8 8 


m STR SCO ; =6.38-T 


( 8 ) 


(9) 


The mass of the concentrator array can be expres$pd as 

m array( T ) -6-34-T 
The dry mass of the vehicle can then be calculated as 

m dry STRl 01 prop> TF .' r > p bus* a ) : = TF m prop -|- m e | ect (a , P bus ), +- m str( t ) + m array (T) ( 10) 

8 8iVen Pr0pella "‘' mass or,he vehide “ 

m b_STR( m prop’ TF . T .I > bus' nl p|.»)->>'dry_STR( lll pn J p. TF .T,P bu!i .o) + m p| (11) 

M!? in , B d he ““ ma 7 Spa “ craft ““ PPoP-kion system involves an itertnive 
procedure. First an educated guess is made of the required mass of the propellant Using this guess the 
initial mass can be calculated by the rocket equation. using tnis guess me 


m o( m prop- TF .T.P bus ,a.m p| .Av.I sp , n ) :=m b _ STR (n, prop ,TF > T,P bus ,m pl ,a).expf_- 

l St -I 


Av 


— rl (12) 


A , r sp^MHD ( T »^ sp bus tT l) i 

th,Tn W t V f Ue f ° r thC m f 5 ° f the P r °P ellant can be obtained by realizing that the mass of the propellant is 
mPtS maSS thC maSS 31 bum0Ut - Mathcad ’ s fi n d function, which uses the Levenberg-MarquJrdt 

tanLcP r aS t USed . C ° nverge on a flnal va,ue for the mass of *0 propellant. This formulation 8 lets us vary 

age fraction, bus power, specific power, payload mass, Av, I sp , and MHD efficiency. 

reqUired f ° r ^ tranSfer iS C3,CUlated by diYiding the required mass of Pliant by the mass 
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t 


Transfer 


(13) 


;= m prop i p bus 

m dot i^J S p ) 

Calculating the transfer time in this manner assumes a continuous bum. 

plane changes are required, such as LEO-GEO transfers with a 28 5° planed change AAvnfZ lI/T § 
used thts mission nude study for LEO - GEO transfers, which 00^1^11^1™ 

LEO to GEO Orbit Transfer 

■nte LEO to GEO orbit transfer mission defined for this study delivers a 5,000-kg payload from low earth 
.hi le £ 32 ^ 7 * 4V ' S “ re • *«*» *• 1st to mat:,, io oft 

was tstatljsela^ot ™ ^ M ’“ ! ™ ise »f *« MHD system 

tkett'fif m MotZr e ?? 0n u affeCtS ,he '” i,ial ma “ of vehide - the I, of the solar thermal 
J™, !u , ? 6 ° f.? 0nds and lts thrust vaned - This is equivalent to varying the mass flow rate of 

he solar thermal rocket. The variation of the initial mass as a function of the bus power or augmented I is 
shown m Figure 2 for unaugmented thrust of 1, 5, and 10 Newtons. As we IncretTe Xu* mm2 few 

which MHn* y decrease 'h e ttt. ass savings. Thus, there exist a maximum mass flow rate or thtust above 
whtch MHD augmentatton wtll mcrease the initial mass required to deliver the payload m iKdesTred orb™ 

K™d ht ^ = 5000 k8 ' ” =20% - « - maximuin thtust 

is around 7.5 Newtons, which corresponds to a mass flow rate of 0.9 gm/sec This maximum occurs 

because the power supply required to obtain a sufficient increase in the I® at Imger ma^fiow mtes?s m 

large for the mass of propellant saved to compensate for its mass as shown in Figure I 
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T = 10 N 

Figure 2: Initial mass as a function of the I sp of the augmented STR and the thrust of the STR only. 


Flow rates, smaller than the maximum possible flow rate, have an associated optimal I_ or electrical bus 
power Figure 2 shows an optimal I sp of 1700 sec for T str = 1 Newton and 1030 sec for T str = 5 Newtons. 
Note that there is no optimal I sp for T str - 10 Newtons, because it is above the maximum thrust or flow 
rate. The decrease in initial mass due to MHD augmentation is 4.2 MT for the 1 Newton STR and 0.45 
Newtons for the 5 Newton STR. The power required to achieve the optimal I sp for the 1 Newton STR is 
found from Figure 4 to be 62 kW, while the augmented thrust is increased to 2 Newtons as shown in 
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Thrust of MHD Augmented Rocket 






— n = 20% 


rj = 50% 

— n = 70% 

Figure 4: Specific Impulse of the MHD augmented STR versus bus power. 



T - ION 

Figure 5: Relationship between initial mass and trip time for several values of thrust. 

Figure 2 shows that a significant mass saving is realized through MHD augmentation; however, is the trip 
time required to transfer the payload to geosynchronous reasonable? Figure 5 shows that the trip time for 
the 1 Newton augmented STR is 320 days at the optimum I sp . While this cuts the transfer time by more 
than a factor of two from that of the unaugmented 1 Newton STR, it is still unreasonable for most missions. 


7 




Increasing the bus power further will further decrease the trip time, however, if shorter trip times are the 

kW 'SSTt the " WC 3re “l Uch better ° ff putting a11 of the P° wer int0 thermal propulsion If the 62 
of electrical power required to reach the optimal I sp were used in a STR, an achievable thrust of 50 

S “ " Ut tbe transfertime t0 under 15 da ys- So the quick truth is that the battle between long 

technology"* maSS S3Vm8S Cann °‘ be S ° lved With MHD au 8 mentation using current power supply 


An improvement in either the efficiency or power supply specific mass will result in reduction in both trio 
time and initial mass in LEO. Figure 6 illustrates this for the efficiency of the MHD accelerator Curves 
are shown for efficiencies of 20, 40 and 70%. Thrust of the STR is varied from a SSnIs 

optimum ° n f 3 ° n8 h eaC ^ CU T\ mit,al m3SS ' S P ' 0tted versus frip time corresponding to the 
vaW th f fr 3Ch Va U t °u f thnJSt AS expected the tnP t'me and initial mass decrease for all thrust 
lues as the efficiency of the MHD accelerator is increased. The effect of the maximum thrust for 
effective augmentation can be seen in the rolling over of the curves for 20% and 40% efficiency The 
maximum thrust for effective augmentation at 70% efficiency is around 26 Newtons; therefore, it does not 
roh over in Figure 6. Finally, if an arbitrary limit of a 100 day trip time is imposed, the efficiency of the 
MHD augmentation must be greater than 20% to be beneficial. Thrust levels greater than 3 Newtons and 5 
Newtons result ,n trip times less than 100 days for efficiencies of 70% and 40% respectively The power 

67 q ?kWf° aC fr eVe f' P l timeS ’ While ° perating at optimal V ^ given in Figure 7 as 62.4 kW and 
67.7 kW for efficiencies of 70% and 40% respectively. 



XXX n = 40% 

©OG r| = 70% 

Figure 6: Initial mass versus trip time for optimum I sp conditions. 
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XXX n = 40% 

©GO r) = 70% 

Figure 7: Bus power versus trip time. 

It is also instructive to evaluate performance for different mission Avs. Figure 8 shows the mass saving that 

can be realized for Av s of 4.2, and 8 km/sec. The mass saving as a function of trip time and thrust is also 
given for a Mar’s mission in Figure 10 . uirusi is aiso 



Av = 6 km/sec Trip Time (days) 

~~ Av = 4.2 km/sec 
Av = 8 km/sec 


Figure 8: Trip time vs. initial mass for Av of 4.2, 6 and 8 km/sec. 





n = 7o% 
XXX r) = 40% 
Q<3>6 n = 20% 


Figure 9: Initial mass versus trip time for augmentation of a solar thermal steam rocket with an I of 300 

sec and a futuristic 5 kg/kW specific power. P 

Augmentation also shows a much greater mass saving when augmenting rockets with a lower initial I 5n 
For an advanced power supply with a specific power of 5 kg/kW, the initial mass can be cut in half while 
the trip time is doubled. Although 5 kg/kW is optimistic it should be noted that a propulsion system 
utilizing water has options that other systems do not. First the tank does not have to be kept thermally 
isolated as in the case of cryogenic propellants. It can be used as a heat sink and its surface can be a 
P^ or ' Unique power supplies that use superconducting magnets to supply the magnetic field for the 
MHD system and store energy may be used for a pulsed system. 
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**X n = 4o% 
ooe n = 70% 


Figure 10: Initial mass versus trip time for a Mar’s mission with a Av of 14.5 km/sec. 


Quasi-one-dimensional Analysis 

The system study showed that to achieve a trip time less than 100 days an efficiency greater than 40% must 
be obtained. A quasi-one-dimensional analysis will give the minimum plasma conductivity needed to 
achieve this efficiency. The conversion efficiency of a MHD accelerator can be defined as (Macheret, 
Miles, and Nelson) as the ratio of the push work to the joule heating in the channel. For a Faraday 
accelerator where jx equals zero this reduces to 


r| n 


;= l( n c) ' u B z 
j y -MTl(n e ) u B 5 


Where u is the flow velocity, Bz is the applied magnetic field and jy is the current in the channel. The 
electrical conductivity as a function of electron number density is can be approximated by 


n e q 


e^e 


m e c e( n n'Q en + n e‘Qel) 
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Where ne is the electron number density, nn is the neutral number density, ce is the electron thermal 
ve oc.ty and me ts the mass of the electron. Qen and Qel are the electron - neutral and electron ion 
collision cross sections. Using these relationships the efficiency of an ideal Faraday channel is shown in 
rigure 1 1 as a function of number density. 



Figure 1 1 : Faraday channel efficiency as a function of electron density for an accelerator operating at 0. 1 


To obtain an efficiency greater than 40% requires electon density over 6*10 A 21 e/m3. In this analysis the 
electron collision cross section of the gas is assumed to be 10 A -15 cm A 2 and the potasium seed 40* 10 A -15 

c ™ 2 ' A Pr m thermodynam ! c e 9 uil,lbrium must be at very low pressures and high temperature to achieve 
this level of conductivity. It is hoped that microwave will create a non-equilibrium plasma that will obtain 
the same electron density at much lower gas temperatures. 


Conclusions 


MHD augmented propulsion offers unique performance characteristics. If non-equilibrium plasmas can be 
created and sustained then over 50% mass saving may be realized. When one considers adavances in both 
permenant magnet field strengths and high temperature superconducting magnets, examining MHD 
augmented propulsion systems becomes attractive. 
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